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Abstract: Based on the simulated annealing, we propose an efficient algorithm for generating a low- 
cost multicast routing with delay constraints. The algorithm starts with a backup-paths-set 
from the source node to each destination node by using the Dijkstra shortest path algo- 
rithm, and then generates the corresponding neighborhood structure. When temperature 
decreasing, a new solution is selected from the neighborhood structures according to the ac- 
ceptance probability and replace the old solution. The simulation shows that our algorithm 
is efficient for some actual networks. 
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1 = Introduction 


As computer networks become faster, new multimedia services such as videoconferencing 
are being offered. As a typical multimedia application, videoconferencing requires the multicast 
routing algorithm!]. A good survey for multicast routing algorithms is provided by Salamal2!. 
Multicast routing algorithms can be classified into two categories. The first category is the 
shortest path algorithms, which minimize the cost of each path from the source node to a 
multicast group member node. The other category is the minimum Steiner tree algorithm 
whose objective is to minimize the total cost of the multicast tree. The first heuristic for delay 
constrained minimum Steiner tree problem was given by Kompella!! and is referred to KPP. 
The heuristic is dominated by computation of a constrained closure graph which takes time 
O(ul|v|?). When the link delays and u takes non-integers, KPP is guaranteed to construct a 
constrained tree if it exists. The bounded shortest multicast algorithm, BSMA, is another delay 
constrained minimum Steiner tree algorithm that is constrained the best in terms of tree cost (4), 
BSMA iteratively replaces the edges in the tree until the tree cost cannot be further reduced. 
The time complexity for BSMA is O(k|v|* lg |u|), where k may be very large in the case of large, 
densely connected networks, and makes it difficult to achieve acceptable running time. 

We use one of the well known artificial intelligence techniques, simulated annealing (SA), for 
solving this problem, which has been successful on similar difficult combinatorial problems®7). 
Basically, the algorithm starts from an initial feasible topology and continues its search around 
the neighbor topologies until a certain stopping criterion is satisfied. Its stochastic properties 
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could prevent it from getting stuck to local minimum unlike the traditional greedy local search 
techniques. In this paper, we propose an efficient algorithm based on SA for generating a low- 
cost multicast tree subject to delay constraints. We call this algorithm simulated annealing 
for delay-constrained low-cost multicast routing algorithm (SADLMA). The algorithm starts 
with a backup-paths-set from the source node to each destination node using the Dijkstra Kth 
shortest path algorithm, and then generates the corresponding neighborhood structure. When 
temperature decreasing, a new solution is selected from the neighborhood structures according 
to the acceptance probability and replace the old solution. The simulation shows that our 
algorithm is efficient for some actual networks. 


2 Problem formulation 


Mathematically, the delay-constrained low-cost multicast routing problem can be formulated 
as follows. Given a graph G = (V, E) with node set V and edge set E, we define two objective 
functions, c(u,v) and d(u, v), on each edge (u,v) € E. Let c(u, v) be the cost of edge (u, v) and 
d(u,v) be its delay. We assume that c(u,v) = c(v,u) and d(u,v) = d(v,u). Let P(u,v) be a 
constrained shortest path from u to v, where a constrained shortest path between nodes u and 
v is defined as the path with the least cost from u to v subject to the constraint that the total 
delay along this path from u to v is less than the delay constraint A. Let C(u,v) and D(u, v) 
be the cost and delay of this path, respectively. Let |V| = n, |E] =m. On this graph, we have 
a source node s, and a set of destination nodes M called the multicast group. The vertices from 
the set V — M — s are called Steiner vertices. We try to construct such a Delay Constrained 
Steiner tree T rooted at s that spans the destination nodes in M so that for each node m in 
M, the delay on the path from s to m is bounded by a delay constraint A. 

For each m € M, if p(s, m) is the path in T from s to m, then 


X0 dle) <A, (4) 
e€p(s,m) 
where A is the delay bound value. 
The minimum cost delay constrained multicast tree is a delay constrained Steiner tree T 


such that 
S ce) (2) 


ecT 
is minimised. 


3 Simulated-annealing-based QOS multicasting algorithm 


3.1 The algorithm details 

SADLMA belongs to a source-based routing algorithm, because it assumes that sufficient 
global information is available to the source. 

Step 1 Encoding and the initial solution 

In SADLMA, a solution is encoded as an array of m = |M| elements, where each element is a 
path from source s to a destination node m; € M in multicast tree T, i.e., £ = (pi, P2, -** , Pm); 
where pi = p(s,mi), m EM, 1 <i<k. 
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The initial solution is constructed by randomly select a path from each backup-paths- 
set (28-11). 

Step 2 Backup-paths-set 

For each destination node m; € M, we compute least-cost paths from s to m by using the 
Dijkstra Kth shortest path algorithm to construct a backup-paths-set for generating neighbors. 
Let p; be the path set for destination node i: 


pi= {pl,.…, pi, , pe}, (3) 


where pi is the jth path for destination node i. 
Step 3 Evaluation function 
The evaluation function is used for evaluating state of search. In our algorithm, we take the 


objective function as an evaluation function 
f(z) = Cost(z). (4) 


Step 4 Neighborhood structure 
Let p; be the path set for destination node i, that is 


pi = {ph e, pi, +>, pk}. (5) 


In the selection process, randomly select pË € pi for the destination node i. 


Then obtain the initial solutions 


To = {p}, , ph, pie }, { +++) jaa fm}, for short (6) 
where 1 < jiye sjitte ,jm < k, and ji,--+ ,Ji,--- ,jm are positive integer. 
Randomly select 7 … ,js;5°** ,js,,, Where 1 全 


In the algorithm, we construct neighborhood structure asl8-11] 


N(To) > {{7s Ie Peis rim}; RER ie: y A jm} ty dist? ,js }}. (7) 


Step 5 Initial temperature and decrementing the temperature 

SA needs to start from a high temperature to, so as to allow acceptance of any new neighbor. 
In SADLMA, set to = 200. We use the simplest and most common temperature decrement rule, 
ti+1 = À- ti, where A is a constant close to 1, but smaller than 1. This exponential cooling 
scheme (ECS) was first proposed by Kirkpatrick et all!2]. In our algorithm, A is set to be 0.95. 

Step 6 Acceptance probability 

Acceptance probability fi;(t) controls accepting a new solution T; from neighbors N(T;) of 
current solution Tj, i.e., 


1, cost(T;) > cost(T;), 
fiz = (8) 
exp(—A,;/t), cost(T;) < cost(T;), 


where t is the current temperature, Aij = cast(T;) ~ cost(T;). 
Step 7 Termination rule 
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We use a fixed iteration number as a stopping critertion, and the maximal iteration number 
is 100. The pseudo code for SADLMA is as follow: 


Procedure SADLMA(G=(V,E), s, M, A, 6, c, d) 
Tbest = Tnow = To; iter = 0; t = to; 
Generate backup-paths-set by using k-SPA; 
While iter < MaxIterNum do 
Metronum = 0; 
While metronum < MaxNum do 
GetNeighbors (N(Tnow), Tnow, S, M, pathset); 
Random select T;, where T; € N(Thow); 
If Ac = cost(T;) — cost(Trow) <O then  Thow = T}; 
If cost(Tnow) < cost(Thest) then 
Tbest = Thow; metronum = 0; 
Else if exp(—Ac/t) >random(0,1) then Thow = Tj; 
12. Metronum 十 十 ; 
13. End of while loop 
14. t = à- t; iter+4; 
15. End of while loop 
16. Return Thest 
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3.2 The complexity analysis 

Theorem 1 The time complexity of SADLMA is O(kmn3 + MoImk), where m is the 
group size, n is the network size, and k is the parameter in k-SPA. 

Proof The time complexity of constructing backup-paths-set by using k-SPA is O(kmn*). 
In the worst case, generating neighbors costs O(mk), and the complexity of whole simulated 
annealing is O(MoImk), where Mo and I are maximal number of iterations for constant tem- 
perature and termination rule, respectively. The time complexity of SADLMA algorithm is 
O(kmn3 + MoImk). 


4 Simulation results and discussion 


The SA algorithm described in this paper has been tested on several randomly generated 
networks based on the Waxman’s algorithm!. 

In the first set of the experiments, SADLMA is compared with BSMA and LD for cost 
performance, where LD is minmum delay Steiner tree without considering the tree cost!?:!0 
and BSMA is proved to be within 0.07 of the optimal for small networks!!! Figure 1 shows the 
tree cost for varying network size with the group size being 4 and 6, respectively; average node 
degree of network is 3; and A = 0.3. The source nodes and destination nodes are varying in 
each experiment. It can be seen from Figure 1 that our algorithm is of a better cost performance 
than the LD algorithms, and is close to BSMA, and could construct low-cost trees which satisfy 
the given delay bound and manage the network resources efficiently. 
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(a) group size=4 (b) group size=6 
Figure 1: Tree cost versus network sìze for A = 0.3; node degree=3 
In the second set of the experiments, Figure 2 shows the tree cost for different network size 


with the group size being 5; average node degree of network being 4 and 5, respectively; and 
A =0.3. In general, SADLMA has good cost performance and is feasible and effective. 
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Figure 2: Tree cost versus network size for A = 0.3; group size=5 


Finally, we consider the iteration time of the algorithm. Figure 3 shows the tree cost for 
varying the iteration times networks with number of network nodes being 30 and 40, respectively 
and group size being 4,5 and 6. As shown in Figure 3, the algorithm converges quickly, and 
has desirable characteristics of approximation iterative heuristics, which satisfies the real-time 


requirement of multimedia network. 
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Figure 3: Tree cost versus iteration times for example network 


5 Conclusion 


A routing algorithm (SADLMA), in which the simulated annealing algorithm is applied 


to support multimedia group communication, is suggested. The simulation shows that our 


algorithm is efficient for some actual networks. The algorithm can guarantee the requirement 


of multimedia group communication for quality of service. 
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基于 模拟 退火 法 的 组 播 路 由 算法 
A WA, HFR? 
(1- 中 国 计 量 学 院 数学 系 ， 杭 州 310018; 2- 大 连理 工大 学 数学 系 ， 大 连 116023) 


搞 要 : 本 文 提出 了 一 种 基于 模拟 退火 算法 的 延 时 约束 最 小 代价 组 播 路 由 算法 (SADLMA)。 首 先 ， 本 算 
法 使 用 Dijkstra 第 K 最 短路 算法 建立 了 从 源 节点 到 每 个 目的 节点 的 候选 集 。 然 后 生成 了 相应 的 邻居 结 
构 。 当 温度 下 降 时 ， 根 据 接收 概率 从 邻居 结构 里 把 新 解 选择 出 来 ， 并 且 代 替 旧 解 。 仿 真 试验 表明 本 算法 
对 实际 网 络 是 有 效 的 。 

关键 词 : 组 播 ， 模拟 退火 法 ， 延 时 约束 


